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The full time-dependent four-wave mixing polarization in quantum dots is microscopically calculated, tak- 
ing into account acoustic phonon-assisted transitions between different exciton states of the dot. It is shown 
that quite different dephasing times of higher exciton states in pancake anisotropic InGaAs quantum dots 
are responsible for the experimentally observed 1 1 1 double-exponential decay in the photon echo signal. 
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1 Introduction Dephasing of the optical polarization caused by the interaction with lattice vibrations 
(phonons) is inevitable in solid state nanostructures and presents a fundamental obstacle for their applica- 
tion in quantum computing. Exciton dephasing in InGaAs quantum dots (QDs) has been recently studied 
by time-integrated four-wave mixing (FWM) measurements [2 J. For the ground state excitonic transi- 
tion in QDs, the temporal dynamics of the measured nonlinear polarization shows an initial rapid decay 
within the first few picoseconds which is followed at later times by a much slower exponential decay of 
the zero-phonon line Excited states, in turn, have a quite different behavior. It has been recently 
shown Q] that after an initial fast decay (similar to that of the ground state) the excited state polarization 
has a double-exponential decay. A related non-exponential long-time decay has been also observed in QD 
molecules EJ. 

The initial decay has been studied theoretically by Vagov et al. (5l|6|, taking into account the interaction 
of a single exciton state of a QD with bulk acoustic phonons, and using the exactly solvable independent 
boson model (BM) |7 J. These results reproduce well the experimental FWM signal at short delay times [8J. 
However, there is no long-time decay in the BM. In order to go beyond one has to consider the phonon- 
assisted coupling between different excitonic states in the QD. Recently we have developed a microscopic 
approach for the linear excitonic polarization in QDs, taking into account both diagonal and nondiagonal 
coupling between different exciton states (9]llOI . In the present work, we extend our theory to the nonlinear 
optical response of a multilevel excitonic system coupled to acoustic phonons and calculate the FWM 
polarization of the QD excited states. 

2 Optical response on a sequence of ultrashort pulses In the heterodyne technique used in [1-4], 
semiconductor QDs are excited by a sequence of ultrashort pulses: 

£{t)^EiS{t-h)+E2d{t-t2) + ..., h<t2<... (1) 

The full Hamiltonian of the optically-driven exciton-phonon system has the form 7Y(t) = H+Hj (t), where 
H is the exciton-phonon Hamiltonian and Hj{t) — —S£{t) — d£*{t) is the exciton-light interaction. 
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Here d = A*" '^^e exciton dipole moment operator, |0} is the QD vacuum state, and /i„ = 

/ dr 5'„(re =r, r/i =r) is the dipole moment of the given single-exciton state \n) (n = 1, 2, . . . ). 
The optical polarization then takes the form 

v{t-MM, . . . ) = {uKt)d{t) u{t)) - {ul^{h)ul^{t2) . . . d{t) . . . UE^hWEAh)) , (2) 

where d{t) = e*^*(ie~'^*. The particular form of the electric field, Eq. 0, allowed us to write the 
full evolution operator = Texp{— i Jq e*^* Hj{t')e~'^^^ di'} as a time-ordered product of operators 
UE{t) due to each individual pulse. The latter can be calculated analytically in any order of d, giving up 
to first order 

Ue{Q) = e'^(^^'''+^'^) = cos \^JiE\ + i^^^^iE d^ + E*d) + . . . , (3) 

where = l/^"!^- heterodyne technique with a double-pulse excitation, the component 

proportional to i?i i?| is essentially filtered out from the full polarization: 

n 2 

Si*S|PFWM(t,r), (4) 



„ , N sin \2^,E^ 



2\iiEi 



sin \^E2\ 



\,iE\ 

where Pfwm(A, t) = r)(i^(0)(i(t)(i^(0)) is the FWM polarization, t > is the delay time between 
the two pulses. 



3 Excitonic multilevel system coupled to acoustic phonons Neglecting biexcitonic effects, the Hamil- 
tonian of several excitonic states in a QD linearly coupled to acoustic phonons takes the form 1 10| 

n q n.m q 

where hun is the bare exciton transition energy, and Mnm (q) the exciton-phonon matrix element. The 
FWM polarization (defined in Sec. 2) can be written as a product of two infinite perturbation series 

v-^ o o V *^ / ifdt'V{t'} ~iJdt"V(t") \ 

PFWM(i,r)=.^|M„nAi™pe^(""^-"'"*)/(n|7I„ve-^ \n) {m\T e o |m) \ , (6) 

71, m \ I 

where arranges all negative-time operators in the inverse order. The external brackets (...) in Eq. (|6} 
mean the finite-temperature expectation value which is taken over the phonon system and thus mixes 
interaction operators from both series. Following 1 10| we use the cumulant expansion: For each pair of 
exciton states (rt, m) we calculate the cumulant K{t, t) defined as 

((n| . . . \n){m\ . . . |m)) = 1 + P^^^ + P^^) + . . . = c^^[k'^^) + if (2) + . . .} (7) 

where K^^^ — P^^\ /C'^) — p(2) _ [p(i)]2/2, and so on. In numerical calculations we proceed up 
to second order in the cumulant thus taking into account both real and virtual phonon-assisted transitions 
between different exciton states | lOJ . 

In the special case of a single level, the two series in Eq. ^ can be combined into one, extending the 
time integration from — r to t, the interaction and time-order operators being redefined, respectively, to 
V{t) = 9{t)V{t) and T = 7i„v [T) for negative (positive) times. As in the linear polarization (and in any 
higher-order nonlinear response), the cumulant expansion ends already in first order, giving exactly 

PFWM(i,r) = (*-)exp|-i^*dt'y'* dt" {TV )V [t" ))^ = '^^p!''/;!'^ ^ , (8) 

where PLit) = e"^"^-* exp{-i /J dt' dt" {TV{t')V{t"))} is the hnear polarization jT). Equation (|8} 
reproduces the result by Vagov et al. 0, but is derived here in a more straightforward way. 
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Fig. 1 Energies of exciton excited states in Fig. 2 Dephasing rates (a) and FWM amplitudes 

anisotropic pancake InGaAs QDs calculated with (b) of the bright excited states at T = 10 K. Bi- 

(full curves) and without (dashed curves) Coulomb exponential fit of the ensemble-averaged photon echo 

interaction. A and B are the two bright excited states is shown by stars. The splitting distribution function 

(thick lines). is shown by a dashed line. 



4 Dephasing of QD excited states In order to simulate the measurable time-integrated FWM, we con- 
sider an ensemble of QDs of different size and shape. Variations of the QD size are mainly responsible for 
a wide distribution (^ 30meV in InGaAs QDs |T1) of the exciton transition frequency. Due to the phase 
prefactor in Eq. (|6}, the transition frequency distribution leads to the photon echo effect; The FWM polar- 
ization from different dots destructively interfere at all times except those around t — t, and the measured 
signal is Pccho(T) = Pfwm{t, t). In a proper simulation, the latter has to be additionally averaged over a 
distribution of the QD excited states splitting, which is much narrower (^ 3 meV 1 1 )) and is mainly due to 
QD shape variations. 

To do this, we assume a pancake shape of a QD (having a smaller size in z-direction) with parabolic 
potentials for electrons and holes, which are additionally anisotropic in the xy-plane. Here we concentrate 
on exciton excited states formed from x- and y-type electron and hole states. Without Coulomb interaction 
the electron-hole pair has two bright (XeXh and yeVh) and two dark states (xei/h and yeXh)- Their energies 
are shown in Fig. 1 (dashed curves) as functions of the anisotropy splitting A = Ey^y^^—Ex^Xh- Switching 
on the Coulomb interaction adds to the exciton Hamiltonian a block-diagonal matrix, the bright and the 
dark states not talking to each other (due to parity), and brings in an additional energy splitting within each 
doublet (Fig. 1, solid curves). At A = bright exciton states A and B are, respectively, exact symmetric 
and antisymmetric combinations of the former states, the lower state A accumulating the whole oscillator 
strength [see Fig. 2 (b)]. 

Quite explicitly, for these excited states the photon echo can be now written as 

Pccho(r) = ^i\ Ca{t) e-^^^^+fi% Cb{t) e-^ri.- + 2/1^^1 Cab{t) e'^^^+^^'>' cos{uja-^b)t. (9) 

Remarkably, the dephasing rates Ta and F b which show up in the FWM polarization, Eq. (|9j, turn out to 
be exactly the same as in the linear polarization. The lower state (A) which is considerably split off has a 
rather weak dephasing at T 10 K [Fig. 2 (a)] which is mainly due to virtual transitions. The upper state, 
in turn, is always close to one of the dark states and has a much stronger dephasing due to real transitions. 
It rapidly changes with A having maxima when the level distance is close to the typical energy (^ 2 meV) 
of the acoustic phonons coupled to the QD. The photon echo amplitudes shown in Fig. 2 (b) are mainly 
determined by the exciton dipole moments /i^, /ig. The functions Ci(r) in Eq. (|9} change rapidly at short 
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Fig. 3 Photon echo signal of anisotropic QD calcu- 
lated for A = 2.9 meV (thin curve) and the ensemble- 
averaged photon echo (thick curve), at long times fitted 
by two exponents (circles). Inset: the same in the short- 
time scale. 

times, while becoming constants in the long-time limit. These constant final values are nothing else than 
the zero-phonon weights which at low temperatures do not differ too much from unity. 

The photon echo signal, Eq. (|9}, shown in Fig. 3 for A = 2.9 me V (thin curves) have oscillations 
up to 200 ps due to the mixed AB -term. The rapid short-time decay (coiTesponding to the broadband 
in the spectrum) is also seen as lowering of both maxima and minima (see the inset in Fig. 3). We fi- 
nally perform the averaging over an ensemble of QDs having different xy-anisotropy degree: Pecho{T) = 
dA PcchoiT; A) A/Aq, where Aq = 3meV is used. The averaged function Pocho(T) has 

oscillations only in the time scale of 27r?i/Ao (Fig. 3, thick curves). At later times it is well approximated 
by a sum of two exponents (Fig. 3, circles) with decay times 20.8 ps and 922 ps. The amplitudes and half 
decay rates of the bi-exponential fit are marked in Fig. 2 by stars. Both amplitudes and rates are very close 
to those of the photon echo signal i-'ccho(''') of the particular QD with A — 2.9meV. Thus, we believe that 
the measured double-exponential long-time decay |T| is due to the two bright excited states in QD, which 
have quite different dephasing times at low temperatures. 
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